Finite- Temperature Properties of Attractive Three-Component Fermionic Atoms in 

Optical Lattices 
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We investigate the finite-temperature properties of attractive three-component (colors) fermionic 
atoms in optical lattices using a self-energy functional approach. As the strength of the attractive 
interaction increases in the low-temperature region, we observe a second-order transition from a 
Fermi liquid to a color superfluid (CSF), where atoms from two of the three colors form Cooper 
pairs. In the strong attractive region, we observe a first-order transition from a CSF to a trionic 
state, where three atoms with different colors form singlet bound states. A crossover between a 
Fermi liquid and a trionic state is observed in the high-temperature region. We present a phase 
diagram covering zero to finite temperatures. We demonstrate that the CSF transition temperature 
is enhanced by the anisotropy of the attractive interaction. 

PACS numbers: 03.75.Mn, 67.85.-d, 71.10.Fd, 03.75.Ss 



The study of cold fermionic atoms has attracted con- 
siderable attention. Fascinating aspects of many-body 
effects have been revealed in various phenomena. A 
crossover between a BCS-type superfluid and Bose- 
Einstein condensation (BEC) was observed by control- 
ling the strength of the attractive interaction using Fes- 
hbach resonances 1, 2, 3, 4, 5, 6]. A superfluid-insulator 
transition was observed for ^Li fermionic atoms in an op- 
tical lattice 01 • These phenomena with a highly tuned 
attractive interaction may be difficult to observe in con- 
densed matter physics. One can expect further that the 
cold fermionic atoms show the novel phenomena beyond 
those observed in the condensed matter. It has been 
shown that three-component (colors) fermionic atoms ex- 
hibit characteristic features. Their properties in optical 
lattices at zero temperature have been studied theoreti- 
cally [1, [13, [llj . It was argued that for atoms with a 
weak attractive interaction two of the three colors form 
Cooper pairs, yielding a color superfluid (CSF). As the 
strength of the attractive interaction increases, there is a 
quantum phase transition from the CSF state to the tri- 
onic state, where three atoms with different colors form 
singlet bound states [13, [iH ■ 

In contrast to the detailed investigations that have 
been undertaken at zero temperature, little information 
is available about the finite-temperature properties. Re- 
cently, fermionic atoms with a balanced population of 
three different hyperfine states were successfully created 
[l^ . The temperature at which these atoms were realized 
was T/Tp ^ 0.37[l3, where Tp is the Fermi temperature. 
Studies of the finite-temperature properties are thus in- 
dispensable. In particular, the stabilities of the CSF and 
trionic phases against thermal fluctuations are important 
in terms of observing these novel states in experiments. 

In this paper, we investigate fermionic atoms with 
three different colors (a = 1,2,3) in optical lattices 
at zero and finite temperatures. Using a self-energy 



functional approach (SEA) [H, [ll], we elucidate char- 
acteristics of the CSF, trionic, and Fermi hquid states, 
and study the phase transition and crossover between 
them. In accordance with the conventional model for cold 
atoms in optical lattices [15| , we set the nearest- neighbor 
hopping and the on-site attractive interaction between 
the atoms with different colors. Low-energy properties 
can be thus described by the attractive three-component 
Hubbard model. 



<i,j> a = l 



(1) 



where c\^{cia) is the fermionic creation (annihilation) 
operator for the state with color a in the ith site and 
nia — c\^Cia. t denotes the nearest- neighbor hopping 
integral, fXa is the chemical potential for the atom with 
color a, and Uap{< 0) is the attractive interaction be- 
tween two atoms with colors a and [3. The condi- 
tion /Iq = (\/2)(Uaj3 -\- Ua-y) is imposed to maintain a 
half-filled system with a balanced population of three- 
component atoms. We assume a homogeneous optical 
lattice and neglect the confinement potential for a first 
approximation. 

The SEA allows us to deal efficiently with zero- and 
finite-temperature properties concerning the phase tran- 
sition driven by correlation effects. For Mott transitions 
in correlated electron systems, precise results have been 
obtained as regard to thermodynamic quantities, exci- 
tation spectra, a nd p hase diagrams including the order 
of the transition [lE [3, [Hi- The SEA is based on the 
Luttinger-Ward variational method. This method en- 
ables us to introduce a proper reference system T^ref, 
which has to include the same interaction term as that 



2 



of the original Hamiltonian ([T]) . A simple reference sys- 
tem is explicitly given by the following Hamiltonian, 

Href = J2i ^rcf ' 
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where a]^ {aa ) is the fermionic creation (annihilation) op- 
erator with color a connecting to the ith site in the orig- 
inal lattice. We then obtain the grand potential as 
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where Vt-cci and Srcf are the grand potential and the self- 
energy for the reference system, and trcf and t are pa- 
rameter matrices of the noninteracting terms of the ref- 
erence Hamiltonial and the original Hamiltonian, respec- 
tively. By choosing the parameter matrix t^of to satisfy 
the condition 9f2/9tref = 0, we obtain a proper refer- 
ence self-energy Sref ) which approximately describes the 
original correlated system. Under the conditions for /Xq, 
mentioned above, we can fix the parameters Cca = and 
= A*Q. We choose a particular gauge for the pair- 
ing field so that only A = A12 = — A21 components of 
the variational parameter Aq/3 are nonzero [8,, JJi] ■ From 
the condition 9ri/9A — 0, we obtain A and can discuss 
whether the color superfluidity appears in the original 
system. To take the trionic state into account, we set 
Vi ^ V2 = = V. If the condition d^/dV = is 
satisfied ai V — the system is driven to the trionic 
phase expected to appear in a large \U\/t. By searching 
the optimized parameters V and A, we discuss a CSF- 
trion phase transition in the attractive three-component 
Hubbard model. In the following, we adopt the density 
of states for an infinite dimensional Bethe lattice, which 
is independent of color. The hopping integral is scaled 
as t/\fd to reach a meaningful limit [13] ■ In this non- 
trivial limit, local correlation effects are known to cause 
various interesting phenomena in correlated electron sys- 
tems, such as a Mott transition. Therefore, the infinite 
dimensional systems pave the way for a detailed study of 
correlation effects. 

We calculate the CSF order parameter $ = (cJi42)j 
the quasiparticle weight Z^, and the entropy per site 
S/L ~ —d{il/L)/dT, where L is the number of lat- 
tice sites. The quasiparticle weight Za is inversely pro- 
portional to the effective mass of the fermionic atom. 
We also calculate the single-particle excitation spectra 
(SPES) pa{iLj/t). We do not consider the possibility of 
the trionic density-wave state induced by spatial fluctu- 
ations. 

We begin our discussions with an isotropic interaction 
system at zero temperature. Figure [l](a) shows the re- 
sults for U12 = U23 = U31 = U . In the isotropic in- 
teraction system, Za becomes independent of a. As 
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FIG. 1: (Color online) (a) \U\/t dependence of the CSF order 
parameter "l", the quasiparticle weight Z, and the entropy 
per site S/L. (b) The single-particle excitation spectra of the 
a = 1 atoms pi (uj/t) [= p2{ijj/t)] and those of the a = 3 atoms 
P3,{ijj/t) at U/t = —1.4. (c) The single-particle excitation 
spectra p{u]/t) [= pi{uj/t) — p2{uj/t) — p3,{uj/t)\ at U/t — 
-3.0. 



\U\/t is increased, the CSF order parameter $ increases 
and then takes its maximum value at U/t ~ — 1.9. At 
U/t = —1.94, $ vanishes discontinuously, suggesting the 
first order transition. The quasiparticle weight Z de- 
creases significantly because of the renormalization ef- 
fect and then drops suddenly to zero at U/t = —1.94. 
The entropy S/L jumps to ~ In 2 at the same The 
residual entropy S/L ~ In 2 in the large \ U\/t region indi- 
cates the formation of the localized trionic state at a site. 
The implications of this phenomenon are discussed later 
in this work. We thus conclude that the CSF-trion quan- 
tum phase transition occurs at U/t = — 1.94(= Uc/t), 
which is of the first order. 

We investigate the CSF-trion transition with respect 
to the SPES. In Fig.fHb) and (c), we show the SPES 
for U/t = -1.4 and -3.0. At U/t = -1.4, the SPES of 
colors 1 and 2 have a spectral gap around the Fermi en- 
ergy a; = 0, which is caused by the superfiuid order. Note 
that piiui/t) = P2{io/t). The SPES of color 3 has a Fermi 
liquid peak aX lo — Q. The incoherent peaks at energies 
distant from the Fermi energy are seen not only in the 
Fermi liquid state but also in the CSF state. The results 
indicate that the attractive interaction induces a renor- 
malized CSF state that is consistent with the decrease 
in the quasiparticle weight Z . In our calculations, lead- 
ing excitation processes are taken account of, resulting in 
the multiple incoherent peaks. In the real system, con- 
tinuous broad incoherent spectra are expected to appear 
instead of the multiple incoherent peaks. On the other 
hand, the SPES of atoms with different colors are the 
same at U/t — —3.0, because the trionic state restores 
SU(3) symmetry broken in the CSF state [l3|. The width 
of the spectral gap Atri is close to 2\U\/t, which corre- 
sponds to the energy required to remove one atom from 
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FIG. 2: (Color online) (a) \U\/t dependence of the CSF order 
parameter <E> and the entropy per site S/L. (b) T/t depen- 
dence of $ and S/ L. 



the trionic bound state. 

In Refs. [13, [m, the CSF-trion transition at T = 
was studied below half filling, using a Gutzwiller varia- 
tional calculation combined with dynamical mean field 
theory. It was argued that the CSF-trion quantum phase 
transition is of the second order in contrast with our re- 
sults. The SPA takes precise account of both the high- 
and low-energy properties of the self-energy. Actually, 
in the present system, we obtain the quasiparticle peaks 
in the CSF, the coherent peak in the Fermi liquid, and 
the incoherent peaks as shown in Fig.[ljb). However, in 
Gutzwiller variational calculations the correlation effects 
in the high-energy incoherent peaks are in principle ne- 
glected, although those in the low-energy coherent peak 
are included. When we evaluate the ground state energy, 
both the low- and high-energy contributions play impor- 
tant roles for obtaining precise results. It is considered 
that the accurate handling of the correlation effects in 
the SFA yields the first order quantum phase transition. 

We next investigate the CSF order parameter and 
the entropy at finite temperatures. In Fig.[2)Ja), $ and 
S/L are shown as functions of \U\/t for temperatures 
T/t = 0, 0.05, 0.08, and 0.12. As the interaction strength 
increases for T/t = 0.05 (0.08), $ increases continuously 
and S/L shows a cusp at \U\/t = 0.96 (1.23). The re- 
sults indicate a second order transition from the Fermi 
liquid to the CSF. As \U\/t increases further, <i> and S/L 
jump to zero and a finite value at \U\/t = 1.83 (1.73), 
respectively, which are evidences of the first order tran- 
sition from the CSF to the trionic state. At T/t ^ 0.12, 
$ remains zero in overall |C/|/^- For S/L, no discontinu- 
ity is observed but the peak around U/t — —1.74. The 
results demonstrate that there is a crossover from the 
Fermi liquid to the trionic state around T/t = 0.10. 

We investigate the temperature dependence of $ and 
S/L with the interaction strength fixed at U/t = 
—0.5, —1.3, and —1.9. The results are shown in Fig.l^^b). 
For U/t = —0.5, the superfluidity vanishes except at 
T/t ^ and the entropy exhibits typical Fermi liquid 
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FIG. 3; (Color online) Phase diagram of the isotropic attrac- 
tive three-component Hubbard model at half filling. 



behavior S/L oc T/t. For U/t = —1.3, $ decreases 
continuously with increasing T/t and finally vanishes at 
T/t = 0.09. We see that S/L at U/t = -1.3 clearly 
changes its slope near the transition point T/t ~ 0.09. 
At T/t < 0.09, S/L accords with the sum of the temper- 
ature dependence of the Fermi liquid and the superfluid, 
while S/L at T/t > 0.09 exhibits a Fermi liquid behavior. 
For U/t = -1.9, $ and S/L show jumps at T/t = 0.04, 
indicating a first order transition from the CSF to the tri- 
onic state. At T/t > 0.04, S/L remains almost constant, 
yielding zero specific heat. This feature is a manifesta- 
tion of the trion gap Atri- 

In the trionic phase, the residual entropy at T = is 
S'/L ~ In 2 as seen in Fig.[TJa). The residual entropy is 
caused by pathological behavior typical in infinite dimen- 
sional systems. The effective hopping integral of the trion 
between neighboring sites can be derived as t^s oc /U"^ 
for a large In an infinite-dimensional system, toff 

becomes irrelevant [13] and localized trions emerge. Ac- 
cordingly, there are two degrees of freedom whether or 
not the trion is found in the site even at T = 0, result- 
ing in the residual entropy In 2. In finite-dimensional 
systems, however, foff is relevant and forms a trionic 
Fermi liquid when T < teff- For higher temperatures of 
tcff < T < Atri , the fully mixed states of the thermally 
excited trions are expected to appear and the entropy 
S/L takes In 2 even in finite-dimensional systems. Our 
results for the trionic state may adequately describe the 
characteristic properties in this temperature region. 

By performing the same calculations, we obtain the 
phase diagram as shown in Fig. [31 We find two types 
of transition and a crossover, namely, the second order 
Fermi liquid-CSF transition, the first order CSF-trion 
transition, and the Fermi liquid- trion crossover. The 
crossover line is determined by the peak of S/L. The 
Fermi liquid-trion crossover line starts around the max- 
imum CSF transition temperature T^gp. This result is 
different from that schematically predicted in Ref. (ll| . 
where the trion transition starts at the critical |C/c|/i of 
the CSF phase along the T = fine. It is expected that 
these phase transitions and crossover can be detected 
by the temperature dependence of the entropy shown 
in Fig.[5]Jb). The trionic state may be detected by the 
change in the Fermi surface, since the coherent peak in 
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FIG. 4; (Color online) CSF transition temperature Tcsf as 
a function of \U\/t with R = 1.0,0.5,0.2, and 0. The broken 
lines are the crossover lines for R = 1.0, 0.5, and 0.2 from left 
to right. 

the trionic state disappears as shown in Fig.[ljc) 

Despite the high tunability of the optical lattice sys- 
tems, it is difficult to cool the fermionic system because 
of the Pauli principle. To observe the CSF state, the 
CSF transition temperature has to be increased by con- 
trollable methods. To this end, we direct our attention 
to the effects of the anisotropy of the attractive inter- 
action. In the experiment, the scattering lengths that 
control the interactions between the three atoms are dif- 
ferent In addition, the respective scattering lengths 
are more tunable than temperature in experiments. 

We set U = Ui2 and U' = Un ^1/23, and introduce 
the anisotropy parameter R= U' /U{< 1). Using the cal- 
culations described above, we obtain Tcsf for R = 0.5, 
0.2, and as shown in Fig. 21 We find that Tcsf in- 
creases with decreasing R, although the total strength 
of the attractive interaction decreases for a given |C/|/t. 
When the anisotropy is introduced, the trion formation 



is suppressed. This tendency extends the CSF phase, 
leading to an increase in the CSF transition tempera- 
ture. In the limit i? — s- 0, the system is equivalent to the 
attractive two-component Hubbard model. Superfiuid- 
ity has already been observed experimentally for systems 
adequately described by the attractive two-component 
Hubbard model Comparing the maximum CSF tran- 
sition temperatures in < R < 1, we find that T™^p/t 
of the attractive three-component Hubbard model is at 
least roughly a half that of the two-component model. 
The findings suggest that the CSF transition tempera- 
ture can be controlled via the scattering lengths, and the 
attractive three-component fermionic atom system in op- 
tical lattices is a candidate of observing the novel CSF 
and the trionic state. 

In recent experiments for balanced three-component 
^Li fermionic atoms, an anomalous increase was observed 
in the three-body loss at a certain magnetic field T^l . Its 
origin has been discussed in connection with the Efimov 
trimer states in the continuous system [H, [H, [13] ■ In the 
low-density region of the lattice system, low-energy prop- 
erties can be well described by the continuous model. It 
is thus interesting to discuss the difference and/or anal- 
ogy between the Efimov trimer in the continuous system 
and the trionic bound state in the lattice system. This 
issue constitutes our future study. 
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